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RacoonWW is the first Monte Carlo generator for e + e~ — ► WW — > 4/ (+7) that includes the electroweak O(a) 
radiative corrections in the double-pole approximation completely. Some numerical results for LEP2 energies are 
discussed, and the predictions for the total W-pair cross section are confronted with LEP2 data. 



1. W-PAIR PRODUCTION AT LEP2 

The investigation of W-pair production at 
LEP2 plays an important role in the verification 
of the Electroweak Standard Model (SM). Apart 
from the direct observation of the triple-gauge- 
boson couplings in c + c~ — * W + W~, the increas- 
ing accuracy in the W-pair-production cross- 
section and W-mass measurements has put this 
process into the row of SM precision tests. 

To account for the high experimental accuracy 
[ |l| on the theoretical side is a great challenge: 

the W bosons have to be treated as resonances 
in the full four-fcrmion processes e + e~ — > 4/, 
and radiative corrections need to be included. 
While several lowest-order predictions are based 
on the full set of Feynman diagrams, only very 
few calculations include radiative corrections^bfe 
yond the level of universal effects (see Refs. [ H, Bl 

and references therein). Fortunately, to match 
the experimental precision for W-pair production 
at LEP2 a full one-loop calculation for the four- 
fermion processes is not needed, and it is suffi- 
cient to take into account only those radiative 
corrections that are enhanced by two resonant W 
bosons. The neglected corrections are of the order 
(a/n)(r-w/M-w), i.e. below 0.5% even if possible 
enhancement factors are taken into account. The 
theoretically clean way to carry out this approxi- 
mation is the expansion about the two resonance 
poles, which is called double-pole approximation 
(DPA) [ |J. A full description of this strategy 
and of different variants used in the literature [ 
B> §t 01 ( some °f them involving further approxi- 
mations) is beyond the scope of this article. We 
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Figure 1. Structure of virtual factorizable correc- 
tions with loop corrections in the blobs 



can only briefly sketch the approach pursued in 
RacoonWW [| § . 



2. RADIATIVE CORRECTIONS WITH 

RacoonWW 

In DPA, 0{a) corrections to e + e~ — ► WW — > 
4/ can be classified into two types: factorizable 
and non-factorizable corrections. We first focus 
on virtual corrections. 

Factorizable corrections are those that corre- 
spond either to W-pair production or to W decay. 
Virtual factorizable corrections are represented 
by the schematic diagram of Fig. [I], in which the 
shaded blobs contain all one-loop corrections to 
the on-shell production and on-shell decay pro- 
cesses, and the open blobs include the corrections 
to the W propagators. For the corresponding 
matrix element M. the application of the DPA 



2 




Figure 2. Typical diagram for virtual non- 
factorizable corrections 



amounts to the replacement 
ka = ^(fcw+'^w-) 

R(M^,M^) 

^ (k 2 w+ -M2)(k^-M*y [ > 

where the originally gauge-dependent numerator 
R(k^ + , fc^,_ ) is replaced by the gauge-indepen- 
dent residue i?(M^,M^), and M 2 = - 
iA/wTw is the location of the poles in the com- 
plex k 2 N± planes. The one-loop corrections to this 

residue can be deduced from the known results 
for the pair production [ [lOfl and the decay [ [til] 

of on-shell W bosons. However, the spin corre& 
tions between the two W decays should be taken 
into account. _ __ 

Non-factorizable corrections [ [12L h3|] comprise 

all those doubly-resonant corrections that are not 
yet contained in the factorizable ones, and in- 
clude, in particular, all diagrams involving par- 
ticle exchange between the subprocesses. Such 
diagrams only lead to doubly-resonant contribu- 
tions if the exchanged particle is a photon with 
energy E 1 < Tw; all other non-factorizable dia- 
grams are negligible in DPA. A typical diagram 
for a vixtual non-factorizable correction is shown 
in Fig. H, where the full blob represents tree-level 
subgraphs. We note that diagrams involving pho- 
ton exchange between the W bosons contribute 
to both factorizable and non-factorizable correc- 
tions; otherwise the splitting into those parts 
would not be gauge-invariant. 

In RacoonWW the virtual corrections are 
treated in DPA, including the full set of fac- 
torizable and non-factorizable 0(a) corrections. 

The real corrections are calculated from full ma- 
trix elements for e + e~ — > 4/7, as described in 
Ref. i.e. the DPA is not used in this part. 

In this way, we avoid potential problems in the 



definition of the DPA for the emission of photons 
with energies E 7 ~ On the other hand, 

this asymmetry in the calculation of virtual and 
real corrections requires particular care concern- 
ing the structure of IR and mass singularities. 
The singularities have the form of a universal 
radiator function convoluted with the respective 
lowest-order matrix element squared |A4o| 2 of the 
non-radiative process. Since the virtual correc- 
tions are calculated in DPA, but the full matrix 
element is used for the real photons, a simple 
summation of virtual and real corrections would 
lead to a mismatch in the singularity structure 
and eventually to totally wrong results. There- 
fore, we extract those singular parts from the 
real photon contribution that exactly match the 
singular parts of the virtual photon contribu- 
tion, then replace in these terms the full |.Mo| 2 
by the one calculated in DPA, and finally add 
this modified part to the virtual corrections. This 
modification is allowed within DPA and leads to a 
proper matching of all IR and mass singularities. 
This treatment has been carried out in two differ- 
ent ways, once following phase-space slicing, once 
using the subtraction formalism of Ref. [ [fq] . 

Beyond 0(a), RACOON WW includes soft- 
photon exponentiation and leading higher-order 
ISR effects in the structure-function approach. 
Using G M as input parameter instead of a(0) in 
the lowest-order cross section, also the leading 
higher-order effects from Act and Ap are taken 
into account. On the other hand, in the relative 
correction factor we use a(0) in order to describe 
the couplings of the real photons correctly. 



3. NUMERICAL RESULTS 

3.1. Predictions from RacoonWW 

A survey of numerical results obtained with 
RacoonWW has already been presented in 
Ref. [ H for LEP2 and linear-collider energies. 
Here we only review the W invariant-mass dis- 
tribution given there, extend the results for the 
total cross section, and add some studies of the 
intrinsic ambiguities in the DPA. 

Figure ^| shows the invariant-mass distribution 
of the dti quark pair for the semi-leptonic channel 
e + e~ — > f M /z + du at y/s — 200 GeV at tree-level 
and with electroweak 0(a) corrections for two 
different recombination cuts, M rcc = 5 GeV and 
25 GeV. The recombination of photons with final- 
state charged fermions is performed as described 
in Ref. [ H|: we first determine the lowest invari- 
ant mass M^f built by the emitted photon and 
a charged final-state fermion. If M 7 f is smaller 
than -Mrec, the photon momentum is added to the 
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Figure 3. Invariant-mass distribution of the du 
pair for e + e~ — > z^^ + du and \/s = 200 GeV 
(taken from Ref. [ ||a) 
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Figure 5. Total WW production cross section at 
LEP2, as given by the LEPEWWG [ @ 
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Figure 4. Relative corrections to the invariant- 
mass distribution of the du pair for e + e~ — > 
j^/u+du and y/s = 200 GeV (taken from 
Ref. [ | a) 



one of the corresponding fermion /. The locations 
of the maxima in the corrected line shapes differ 
by up to 30MeV between the two values of M rec . 
As expected, there is a tendency to shift the max- 
ima to larger invariant masses if more and more 
photons are recombined. In Fig. |i| we display 
the relative corrections S — da/da — 1 for the 
two values of M rcc , which illustrates the strong 
dependence of the corrected invariant-mass dis- 
tributions on the treatment of the real photons. 
We obtain consistent results for the phase-space 
"slicing" and the "subtraction" methods. The 
size of the shown effects demonstrates that a 
careful treatment of real photons is mandatory 
in the W-mass reconstruction at LEP2 accuracy. 

Figure H shows a comparison of RacoonWW 
results ana of other predictions with recent LEP2 
data, as given by the LEP Electroweak Working 
Group [ [lj |l6| . The data are in good agreement 
with the predictions of RacoonWW and YF- 
S WW3 [ pf . The predictions of these two gener- 
ators differ between 0.5-0.7%.! More details on 
the conceptual differences of tne two generators, 



t Meanwhile the dominant source of this difference has 
been found, and the new results of YFSWW3 agree within 
0.3% with the results of RACOONWW. Details on the new 
YFSWW3 predictions can be found in Ref. [ | . 
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as well as a detailed comparison of numerical re- 
sults, can be found in Ref. [|3j. Figure || also 
includes the prediction provided by GENTLE [ 
|l~7|l , which differs from the RacoonWW and YF- 
SWW3 results by 2-2.5%. This difference is due 
to the neglect of non-leading, non-universal O(a) 
corrections in GENTLE. Consequently, the com- 
parison between SM predictions with the precise 
measurements of the W-pair production cross sec- 
tion at LEP2 reveals evidence of non-leading elec- 
troweak radiative corrections beyond the level of 
universal effects. 

3.2. Intrinsic ambiguities of the DPA 

In order to investigate the accuracy of the DPA 
quantitatively, we have performed a number of 
tests. The implementation of the DPA has been 
modified within the formal level of qTw/A/wj 
and the obtained results have been compared. 
Because in RacoonWW only the virtual cor- 
rections are treated in DPA, while real photon 
emission is based on the full e + e~ — > 4/7 matrix 
element with the exact five-particle phase space, 
only the 2^4 part is affected by the following 
modifications. We consider three types of uncer- 
tainties (see Ref. [|9| for more details): 

• Different on-shell projections: 

For the DPA one has to specify a projection 
of the physical momenta to a set of mo- 
menta for on-shell W-pair production and 
decay^], in order to define a DPA. This can 
be done in an obvious way by fixing the di- 
rection of one of the W bosons and of one of 
the final-state fermions originating from ei- 
ther W boson in the CM frame of the incom- 
ing e + e~ pair. The default in RacoonWW 
is to fix the directions of the momenta of the 
fermions (not of the anti-fermions) result- 
ing from the W + and W~ decays ("def"). 
A different projection is obtained by fixing 
the direction of the anti-fermion from the 
W + decay ("proj") instead of the fermion 
direction. 

• Treatment of soft photons: 

As explained above, the matching of IR 
and mass singularities between virtual and 
real corrections requires a redistribution of 
the singular parts. This redistribution fixes 
only the universal, singular terms, while 
the redistribution of non-singular terms are 
mere convention. Owing to the asymmet- 
ric treatment of the corrections (virtual in 



§This option only illustrates the effect of different on-shell 
projections in the four-particle phase space. If real pho- 
tonic corrections were treated in DPA the impact of dif- 
ferent projections could be larger. 



DPA, real from full matrix elements), dif- 
ferent redistributions of non-singular con- 
tributions change the result by terms of the 
order (a/7r)(r\y/Mw), i- e - these redistri- 
butions are equivalent within the accuracy 
of the DPA. In RacoonWW two different 
schemes for this redistribution are imple- 
mented. As default, the endpoint contri- 
butions of the subtraction functions, as de- 
fined in Ref. [ , are calculated in DPA 
and added to the virtual photon contribu- 
tion. In the other scheme, the universal 
IR-sensitive part is extracted from the vir- 
tual photon contribution a la YFS [ |T^] and 
added to the real photon contribution. The 
resulting soft+virtual part of the photonic 
correction is, thus, treated off shell ("eik"). 

The difference between the two described 
treatments is that certain terms of the form 
(a/w) x 7r 2 x O(l) are either multiplied with 
the DPA ("def") or with the full off-shell 
Born cross sections ( "eik" ) . 

• On-shell vs. off-shell Coulomb singularity: 
The Coulomb singularity is (up to higher 
orders) fully contained in the virtual 0(a) 
correction in DPA. Performing the DPA 
to the full virtual correction leads to the 
on-shell Coulomb singularity, which is a 

simple factor of cot/ (2/3), where (3 is the 
velocity of an on-shell W boson. How- 
ever, since the Coulomb singularity is an 
important correction in the LEP2 energy 
range and is also known beyond DPA [ 
Q, RacoonWW includes this extra off- 
shell Coulomb correction factor as default. 
This replacement of the Coulomb singular- 
ity is performed by adding and subtract- 
ing the corresponding contributions in the 
virtual non-factorizable corrections, as de- 
scribed in Ref. [ [l3| . Switching the extra 
off-shell Coulomb terms off ("Coul"), yields 
an effect of the order of the uncertainty of 
the DPA. 

In the following table and figures the total 
cross section and two distributions for e + e~ — > 
ud fj,~v p (jy) have been compared for the different 

versions of the^ DPA defined above. The results 
have been obtained using a photon-recombination 
procedure that is similar to the one decribed in 
the previous section. The precise definition can 
be found in Ref. [ |]. 

The results for the total cross section are shown 
in Table 111. We find relative differences at the 
level of fM%. As expected, the prediction that 
is based on the on-shell Coulomb correction is 
somewhat higher than the exact off-shell treat- 
ment, since off-shell effects screen the positive 
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y/s/ GeV 








172 


-0.03 


-0.09 


0.79 


200 


-0.03 


-0.01 


0.13 


500 


-0.01 


0.08 


0.01 



Table 1 

Intrinsic DPA ambiguities 5 = a /add — 1 of the 
RacoonWW predictions for the total cross sec- 
tion of e + c~ — > \idfi~ P^j) (based on the results 
ofRef. [§) 



Coulomb singularity. Note that for the low CM 
energy of 172 GeV the difference between on-shell 
and off-shell Coulomb singularity, which is quite 
large (0.79%), cannot be viewed as a measure 
of the theoretical uncertainty, since the on-shcll 
Coulomb singularity is not adequate near thresh- 
old. 

In Figs. H and U\ we show the differences of the 
"proj", "eiR", and "Coul" modifications to the 
default version of the DPA for the distributions 
in the W-production angle 6\v+ and in the ud 
invariant-mass M u g, both at y/s = 200 GeV. For 
the # w + distribution the relative differences are of 
the order of 0.1-0.2% for all angles, which is of the 
expected order for the intrinsic DPA uncertainty. 
For the M u g distribution the DPA uncertainties 

are at the level of 0.1-0.3% within a window of 
2r-w around the W resonance. The uncertainties 
grow with the distance from the resonance point, 
as expected, since the DPA runs out of its range 
of validity away from the resonance region. 
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